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The  purpose  of  Chis  report  is  to  present  end  describefmethods 
which  may  be  used  to  define  rocket  exhaust  flow  fields,  launch  tube 
internal  pressure  profiles,  recoil  forces,  and  effects  on  missile 
launch  due  to  launcher-missile  interactions. 

To  this  end,  computer  programs  have  been  written  which  will 
calculate  relevant  Bata  for  input  into  other  programs  giving  plots 
of  the  stagnation  pressure  and  temperature  profiles  in  a  given  rocket 
exhaust  field.  A  program  hag  been  (.written  which  will  give  a  plot  of 
internal  tube  pressure  at  the  rocket  nozzle  exit  plane  as  a  function 
of  missile  ^tjravel.  A  method  for  estimating  recoil  is  explained  and 
a  program, has- been  developed  which  simulates  the  launcher  in  the 
pitch  plane  for  use  in  determining  the  launcher  and  missile  motions 
in  this  plane  and  the  effects  on  the  missile  of  thrust  misalignment  and 
center  of  gravity  offsets. 

4 


ii 


FOREWORD 

Acknowledgement  is  made  to  Mr.  Dean  Christensen  for  his 
major  contributions  in  writing  the  subject  computer  programs 
and  to  Mr.  Billy  Campbell  for  his  contributions  in  the  same 
vein. 


ill 


i 


/ 


CONTENTS 

Page 


Section  I.  INTRODUCTION  1 

Section  II.  EXHAUST  FLOW  FIELD  DEFINITION 

2.1  Introduction  .  . .  2 

2.2  Symbol  Definitions . .  .  . .  2 

2.3  Exit  Mach  Number  and  Intercept  of  Linear 

Spreading  Characteristic  Curve  .  .  .  . .  3 

2.4  Axial  Temperature  Decay  Curve  (with  Plot)  ......  6 

2.5  Axial  Pressure  Decay  Curve  (with  Plot)  .  .  10 

2.6  Boundary  Coefficients  (Region  I)  .  .14 

2.7  Nozzle  Placement  .  .  :  .  . . 18 

2.8  Boundary  Plot  (Region  I) .  20 

2.9  Pressure  Plume  Plot  (Region  II)  .  22 

2.10  Pressure  Plume  Plot  (Region  III) . 27 

2.11  Calculation  and  Plot  of  Isotherms  in  Exhaust 

Field  (Region  II) .  ............  32 

2.12  Calculation  and  Plot  of  Isotherms  in  Exhaust 

Field  (Region  III) . . .  39 

Section  III.  TUBE  PRESSURE  PROFILE  AND  RECOIL  PREDICTION 

3.1  Tube  Pressure  as  a  Function  of  Length  ........  44 

3.2  Recoil  Prediction  .  50 

Section  IV.  LAUNCHER  SIMULATION 

4.1  Introduction  . . 52 

4.2  Symbol  Definitions  .  .  .  ...........  52 

4.3  Launcher  Simulation  Equations  .....  .  .  .  59 

4.4  Description  and  Program  ....'.* . .  .  .  .  67 

Section  V.  CONCLUSIONS  113 

REFERENCES  ...........................  114 


% 


V 


Ln 


ILLUSTRATIONS 

Figure  Page 

2-1  Nozzle  Placement,  Region  I  Pressure  Plume 

and  Region  II  Pressure  Plume  .  .....  26 

2-2  Pressure  Plume  Plot  Region  III . 31 

2-3  Isotherm  Profile  Plot  Region  II . 38 

2- 4  Isotherm  Profile  Plot  Region  III . 43 

3- 1  Internal  Tube  Pressure  at  Nozzle  Exit 

Versus  Missile  Travel  .  49 

3- 2  Estimation  Method  for  Recoil  Forces . 51 

4- 1  Launch  Simulation  Missile  and  Launcher  Model  .  54 

4-2  Firing  Position  Models  for  Launch  Simulation  .  55 

Missile  Pitch  Angle  Versus  Time . .  .  56 

Missile  Pitch ‘Rate  Versus  Time  . .  57 

Missile  Y-Velocity  Versus  Time  .  ....  58 

Table 

4-1  Launch  Simulation  Memory  Map  ......  i  .  68 


vi 


Section  I.  INTRODUCTION 


The  object  of  this  task  is  to  provide  parameters  and  information 
relevant  to  the  reduction  in  weight  of  shoulder  fired  rocket  launchers 
to  a  minimum  value  consistent  with  the  maintenance  of  required  system 
performance. 

The  scope  of  this  effort  encompasses  the  following:  (1)  Prepara¬ 
tion  of  computer  programs  to  define  the  stagnation  pressure  and  tempera¬ 
ture  profiles  in  the  exhaust  field  of  a  rocket,  (2)  Develop  a  method 
for  predicting  the  internal  pressure  in  a  launch  tube  during  launch, 

(3)  Discussion  of  means  for  estimating  recoil  forces,  and  (4)  Establish 
a  computer  simulation  program  for  predicting  launcher  effects  on  aiming 
and  launching  of  the  missile. 


Section  II.  EXHAUST  FLOW  FIELD  DEFINITION 


2.1  Introduction 


The  overall  purpose  of  the  rocket  exhaust  flow  field  program  is 
to  provide,  given  all  necessary  data  for  a  rocket,  a  visual  representa¬ 
tion  of  the  stagnation  pressure  and  stagnation  temperature  profiles 
in  the  rocket  exhaust.  Several  programs  are  provided  which  calculate 
input  data  needed  by  the  temperature  and  pressure  programs. 

The  equations  and  methods  upon  which  these  programs  are  based 
are  to  be  found  in  the  references  listed  in  the  reference  section.  Each 
program  will  be  shown  listed  on  programming  forms  with  input  data 
shown  in  the  appropriate  registers  and  output  values  blocked  in  with 
dashed  lines.  A  description  will  be  given  with  each  program  giving 
its  purpose,  inputs,  outputs  and  equations. 

The  machine  for  which  these  programs  are  written  is  the  Hewlett- 
Packard  9100B  with  associated  Plotter  and  Extended  Memory  units. 

2-2  Symbol  Definitions 


b 

c 

De 

Dt 

Ir 

It 

K 

Is 

Mb 

Me 

mp 

mt 

Pc 

q 

qe 

%1S 


^rns 


max 


Region  I  Boundary  Coefficient 
Region  I  Boundary  Coefficient 
Region  I  Boundary  Coefficient 
Nozzle  Exit  Diameter  (inches) 

Nozzle  Throat  Diameter  (inches) 

Intercept  of  Linear  Spreading  Characteristic  Curve. 
Intercept  of  Axial  Temperature  Decay  Curve 
Ratio  of  specific  heats  of  exhaust  products 
Length  of  Supersonic  Cone  in  Exhaust  Flow  (inches) 
Maximum  Mach  Number  at  End  of  Region  I 
Exit  Mach  Number 

Slope  of  Axial  Pressure  Decay  Curve 
Slope  of  Axial  Temperature  Decay  Curve 
Nominal  Design  Chamber  Pressure  (PSIA) 

Value  of  Isobar  to  be  Plotted 

Gage  Stagnation  Pressure  at  Nozzle  Exit  (PSIG) 

Gage  stagnation  pressure  on  axis  at  end  of  supersonic 
cone  (PSIG) 

Maximum  value  of  pressure  parameter  for  plot  (PSIG) 

Exhaust  Plume  Radius  at  End  of  Region  I 
Plot  Scale  Factor  for  X-Axis 
Plot  Scale  Factor  for  Y  -  Axis 
Value  of  isotherm  to  be  plotted  (°F) 

Ambient  Temperature  (°F) 

Flame  Temperature  of  Propellant  (°F) 
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Tmax  “  Maximum  value  of  temperature  parameter  for  plot  (°F) 

Xmax  "  Maximum  value  of  distance  parameter  for  plot 

dX  Plot  increment  distance 

0^  -  Nozzle  expansion  half  angle  (radians) 

2.3  Exit  Mach  Number  and  Intercept  of  the  Linear  Spreading  Characteristic 
Curve 

2.3.1  Purpose 

To  calculate  the  Mach  Number  of  the  Exhaust  Flow  at  the  nozzle 
exit  plane,  and  to  determine  the  value  of  the  linear  spreading 
characteristic  curve  intercept.  This  last  parameter  is  necessary  in 
the  generation  of  the  Exhaust  Pressure  Plume  Plot. 

2.3.2  Inputs 

De ,  Dt ,  K 

2.3.3  Outputs 
Me,  Ir 


2.3.4  Equations 


At  Dt2  r 

—  -  -  Me  — 


2  +  (K-l)  Me2 


n 20 


K+l 

(X-l) 


This  equation  is  solved  iteratively  for  Me- 

Ir  *  0.15489201  Me3-  2.2941041  Me2  +  i0. 612059  Mgi-4. 0070141  '  (Ref  3) 

2.3.5  Description 

After  program  is  entered  the  values  of  K,  Dt,  De  are  entered 
in  the  Z,  Y,  X  registers  respectively.  Press  Continue  and  the  value 
for  Me  will  be  displaytd  in  the  X,  Y,  Z  registers.  Press  Continue 
twice  and  the  value  for  Ir  will  be  displayed  in  the  X,  Y*  Z  registers. 


V 


2.4  Axial  Temperature  Decay  Curve  (with  Plot) 

2.4.1  Purpo  e 

To  generate  values  for  the  slope  and  intercept  of  the  Axial 
Temperature  Decay  Curve. 

2.4.2  Inputs 

Met  De j  Tc ,  Xmax »  Imax 

The  values  for  and  Tmax  should  be  for  a  6  in.  by  9  in.  plot. 

2.4.3  Outputs 

Mt.  It 

Plot  of  stagnation  temperature  on  axis  vs  distance  downstream 
on  the  axis. 

2.4.4  Equations 

_mt  =  .020042599Me3  -  .24382049Me2  +  1.0618553M3  +  .11169953  (Ref  3). 
It  =  6.36  Me  -  2.874  (Ref  3). 

2.4.5  Description 

After  program  has  been  entered,  the  values  of  Mg,  De,  Tc  are 
entered  in  the  Z,  Y,  X  registers  respectively.  Press  Continue  and.  the 
value  for  Mt  will  be  displayed  in  the  Y  register,  value  of  It  in  the 
X  register.  Press  Continue,  then  enter  values  for  Tmax  and  Xma*  in 
Y  and  X  registers  respectively  for  a  6  in.  by  9  in.  Plot.  Press 
Continue  for  Plot. 


To  generate  values  for  the  slope  of  the  Axial  Pressure  Decay 
Curve  and  the  length  of  the  supersonic  cone  in  the  exhaust  flow. 


2.5.2  Inputs 


Me »  De »  Pc  K ,  Xroax ,  9msmax 

The  dimensions  of  De  should  be  in  feet. 


Values  of  X^  and  qms 


should  be  for  a  6  in.  by  9  in.  Plot. 


2*5.3  Outputs 


is.  mp  .  9  e .  9ms 


2.5.4  Equations 


Is,  =  16.113Me  -  16.125  (Ref  3) 
«p  -  -.451667Me  -  1.375  fRef  3) 

■  14.7  | ^  -  lj  (Rel 

fr-  *1*-  r* 

-  *  jpwj  L“ 


2.5.5  Description 


_2K_  2 

K+l  ^ 


■  C) 


14.7  (Ref  3) 


After  program  is  entered,  the  values  of  Mg,  De,  Pc  are  entered 
in  the  Z,  Y,  X  registers  respectively.  Press  Continue,  the  value  for 
1 s  in  feet  will  be  displayed  in  the  Z  register,  *s  in  inches  in  the  Y 
register,  mp  in  the  X  register.  Press  Continue,  enter  value  of  K  in 
X  register,  press  Continue.  The  value  for  9ms  will  be  displayed  in 
the  Z  register ,(  Ijj,/ 9  e)s  in  the  Y,  9e  in  the  X.  Press  Continue, 


then  enter  9ms_ 


in  Y  register,  Xmax  in  X  register  for  a  six  inch 


by  nine  inch  Plot.  Press  Continue  for  Plot. 


Stop 


2.6  Boundary  Coefficients  (Region  I) 

2.6.1  Purpose 


To  generate  coefficient  values  which  determine  the  extent  of  the 
supersonic  exhaust  flow  region  (Region  I). 


2.6.2  Inputs 

Me i  De ,  Dt  >  Pc,  ®e 

2.6.3  Outputs 

a,  b,  c,  Mb,  ri 

2.6.4  Equations 


r^  -  re  sec  w 


1  -  sec  w 


(Ref  1) 


b  *  Tan  w  (re  *  a) 


(Ref  1) 


Qre  -  a)2  +  b2]  7  (Ref  I) 


[l*l¥ 


(Ref  1) 


Cl  -  2  [1+  *t-  (Ref  1) 

2.6.5  Description 

After  program  is  entered,  values  for  Mg,  De,  Dt  are  entered  in 
Z,  Y,  X  registers  respectively.  Press  Continue,  then  enter  Pc,  K,  0e 
in  Z,  Y,  X  registers  respectively.  Press  Continue,  and  values  for  a, 
b,  c  will  be  displayed  in  the  Z,  Y,  X  registers  respectively. 

Press  Continue,  values  for  fy,  rj^  will  appear  in  Y,  X  registers 
respectively. 


Code 


2 . 7  Nozzle  Placement 


2.7.1  Purpose 

To  generate  a  representation  of  a  nozzle  as  a  visual  aid  when 
plotting  pressure  fields  in  the  exhaust. 

2.7.2  Inputs 

®e *  ®e >  ®ti  SFX,  SFy 

9e  should  be  in  degrees. 

2.7.3  Outputs 

Plot  of  nozzle 

2.7.4  Description 

After  program  is  entered,  enter  9e .  De>  Dt»  In  Z,  Y,  X  registers. 
Press  Continue,  enter  SFy,  SFX  in  Y,  X  registers.  Press  Continue 
to  generate  plot  of  nozzle. 


J^Wv'V 


To  generate  a  plot  showing  the  extent  of  Region  I  in  the  rocket 
exhaust  flow  field. 

2.8.2  Inputs 

a.  b,  c,  Ax,  SFX,  SFV 

y  b 

A  good  value  for  ax  is  -Jq 

2.8.3  Outputs 

Plot  showing  the  exhaust  plume  boundary  in  the  supersonic 
region  (Region  I). 

2.8.4  Equations 

C2  =  (r-a)2  +  vx-b)2  (Ref  1) 


2-8.5  Description 

After  program  is  entered,  values  for  a,  b,  c  are  entered  in  the 
Z,  Y,  X  registers.  Press  Continue,  then  enter  Ax,  SFy,  SFX  in  the 
Z,  Y,  X  registers.  Press  Continue  and  a  plot  of  the  Region  I  Boundary 
will  be  generated. 


2.9  Pressure  Plume  Plot  (Region  II) 

2.9.1  Purpose 

To  generate  the,  isobar  plots  in  the  mixing  region  (Region  II). 

2.9.2  Inputs 

Ir,  b,  De,  lg,  r^,  Ax,  SFx,  SFy,  jm^J  ,  9mg,  1  i 

2.9.3  Outputs 

Isobar  plots  for  Region  II  of  the  exhaust  flow. 

2.9.4  Equations 

Ax-  (  18  -  b)  .02  +  1  (10"6) 


1 


2.9.5  Description 


After  program  is  entered,  enter  Ir,  b,  De  in  Z,  Y,  X  registers. 
Press  Continue,  enter  ls,  rj_.  Ax  in  Z,  Y,  X  registers.  Press 
Continue,  enter  SFy,  3FX,  in  Y,  X  registers.  Press  Continue,  enter 
JmpJ,  9ms,  9i  in  Z,  Y,  X  registers.  Press  Continue  to  generate  plot 
of  9^.  After  each  isobar  is  plotted,  enter  value  for  next  9  ^  in  X 
register  and  press  Continue. 
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Distance  Downstream  fr 


2. 10  Pressure  Flume  Plot  (Region  III) 

2.10.1  Purpose 

To  generate  a  plot  of  isobars  in  the  subsonic  exhaust  flow 
region  (Region  III). 

.10.2  Inputs 

^r>  >  De,  mp ,  ls  ,  Xmax>  rmax,  <l»  %is 

Values  of  Xtnax  and  rmdx  should  be  for  a  five  inch,  by  fifteen 
inch  plot. 

2.10.3  Outputs 

Isobar  plots  for  downstream  distances  versus  radial  distances 
from  center  of  nozzle  exit. 


After  program  is  entered,  Ir,  ‘Ig,  0e  are  entered  in  Z,  Y,  X 
registers  respectively.  Press  Continue,  then  enter  mp,  ls,  qms  in 
Z,  Y,  X  registers  respectively.  Press  Continue,  then  enter  values  of 
Xj,,^  and  rmax  for  a  five  inch  by  fifteen  inch  plot.  Press  Continue, 
then  enter  in  X  register  value  of  first  isobar  to  be  plotted.  Press 
Continue  for  plot.  Upon  completion  of  plot,  enter  next  isobar  value 
and  press  Continue  for  plot. 


2.11  Calculation  and  Plot  of  Isotherms  in  Exhaust  Field1  (Region  II) 

2.11.1  Purpose 

To  calculate  and  plot  isotherms  in  the  mixing  region  from  the 
supersonic  cone  outwards. 

2.11.2  Inputs 

Dg  f  Ij- ,  lg,  I  £  g  m  r  i  Tg  i  SFX,  SFy ,  Tg }  ax,  X,  T 

2.11.3  Outputs 

Txs  *  stagnation  temperature  on  axis  at  end  of  supersonic  cone 
Isotherms  will  be  plotted  on  a  radial  versus  axial  distance 

scale. 

2.11.4  Equations 


2.11.5  Description 


This  program  is  composed  of  five  subprograms,  Po 


which  are  er.te 
be  on  with  the 
Inputs  as  shown 
for  T$g.  Enter 
superspnic  coni: 
the  value  of  t  v 


less  than  T. 


xs  ■ 


ed  into  the  Extended  Memory  Unit.  The  Printer  should 
Y  button  depressed.  After  program  is  entered  make  all 
for  P0.  Press  Continue  and  a  value  will  be  printed 
this  value  as  shown  for  ?2  and  the  outline  of  the 
will  be  plotted.  After  each  isotherm  is  plotted  enter 
ie  next  isotherm  to  be  plotted.  These  values  must  be 
The  curves  will  be  drawn  for  the  positive  radius  only. 


Vf. 


Code 


Step 


Code 


m  o 

r-»  ia 


(saqoui) 


CM 
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2.12.1 

Purpose 

in  the 

To  calculate  and 
region  downstream 

plot  isotherms  in 
of  the  end  of  the 

a  rocket  exhaust  flow  field 
supersonic  cone. 

2.12.2 

Inputs  . 

Dei  Ir »  1 s  >  It  > 

mt »  Tc,  SFX ,  SFy , 

Ax,  X,  Ta,  T 

2.12.3 

Outputs 

scale. 

Isotherms  will  be  plotted  on  a  radial  versus  axial  distance 

2.12.4 

Equations 

2.12,5  Description 

After  program  is  entered,  all  necessary  inputs  are  made  as 
indicated  on  program  forms.  For  T,  enter  value  of  isotherm  desired. 
After  it  is  plotted  enter  next  value  of  T  and  continu  this  procedure 
until  plot  is  completed. 
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Step 


Section  III.  TUBE  PRESSURE  PROFILE  AND  RECOIL  PREDICTION 


3. 1  Tube  Pressure  as  a  Function  of  Length 

3.1.1  Discussion  of  Program 

A  method  was  developed  whereby  the  pressure  inside  a  launch 
tube  due  to  missile  exhaust  gases  could  be  predicted  as  a  function  of 
missile  travel.  The  procedure  used  is  discussed  in  the  following. 

A  close  fit  was  assumed  between  the  missile  motor  nozzle  and 
the  walls  of  the  launch  tube,  with  no  space  available  to  allow  the 
exhaust  gases  to  expand  around  the  nozzle  rim  in  a  forward  direction. 

The  pressure  on  the  launch  tube  walls  was  thus  assumed  to  be  due  to 
the  component  of  exhaust  gas  flow  parallel  to  the  nozzle  wall  exerting 
a  force  on  the  launch  tube  walls.  This  program  gives  the  pressure  on 
the  launch  tube  at  the  location  of  the  rocket  nozzle  exit  plane  at  the 
time  the  exit  plane  passes  each  given  location. 

This  procedure  does  not  account  for  shock  wave  buildup  in  the 
tube  behind  the  nozzle  nor  for  any  effects  attributable  to  them. 

This  would  be  a  consideration  for  future  effort  combining  theoretical 
and  test  results. 

3.1.2  Results  and  Discussion  of  Accuracy 

The  program  as  written  will  provide  a  prediction  of  tube  pressure 
versus  missile  travel.  The  accuracy  of  the  program  will  not  be 
ascertainable  until  comparisons  with  actual  test  results  can  be  perform¬ 
ed. 

The  predictions  may  give  lower  than  actual  values  since  the 
effects  of  shock  waves  in  the  tube  were  not  considered.  Shock  waves 
could  cause  pressure  variations  aft  of  the  missile  which  may  exceed 
predicted  values.  If  test  results  seem  to  indicate  this,  a  study 
could  be  undertaken  to  simulate  shock  wave  effects. 

3.1.3  Tube  Pressure  at  Nozzle  Exit  Plane  Versus  Missile  Travel 

3.1.3. 1  Purpose 

To  calculate  and  plot  values  for  the  pressure  insidje  a  launch 
tube,  due  to  exhaust  gases,  at  the  nozzle  exit  plane  as  a  function  of 
missile  travel. 

3. 1.3.2  Inputs 

W  -  Effective  Missile  Weight  (Pounds) 

Pr  "  Pcl/thrus^ 
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1 


Pxy 

-  (  +  14. 7)/Pc 

©n 

Nozzle  half  angle 

cn 

Time 

Tn 

Thrust 

At 

Time  increment 

SF„ 

Plot  ordinate  scale 

factor 

Plot  abscissa  scale 

factor 

3. 1.3.3  Outputs 

Plot  of  tube  pressure  (PSIG)  versus  missile  travel  in  tube,  (inches). 


3. 1.3.4 


After , program  is  entered,  enter  PXy,  Pr,  W  in  Z,  Y,  X 
registers.  Press  Continue,  enter  9n,  SFy,  SFX,  in  Z,  Y,  X  registers. 
Press  Continue,  enter  <k-  in  X  register.  Press  Continue,  enter  Tn,  tn, 
in  Y,  X  registers.  Press  Continue,  value  of  (tube  pressure,  missile 
travel)  will  be  plotted,  continue  entering  values  for  Tn,  tn  until 
plot  is  completed. 


Step 


Step 


Figure  3-1  .  Internal  Tube  Pressure  at  Nozzle  Exit  Versus  Missile  Travel 


3 . 2  Recoil  Prediction 


The  recoil  force  predictions  are  based  on  the  following 
considerations:  missile  shoe  frictional  forces,  gas  friction  forces 
on  the  walls  of  the  launch  tube,  gas  impingement  on  the  internal 
step  for  a  stepped  telescoping  launcher  and  gas  impingement  on  the 
forward  face  of  the  tube  at  missile  exit.  These  forces  all  contribute 
to  a  total  force  aft  (recoil)  during  a  missile  firing.  Any  internal 
protrusion  other  than  a  step  would  also  have  to  be  accounted  for  in  a 
recoil  estimate. 

The  launch  simulation  program  described  in  Section  IV  can  be 
used  to  obtain  the  force- time  histories  of  the  shoe  friction  forces 
and  the  time  at  which  the  exhaust  gas  would  begin  to  impinge  on  the 
step  or  any  other  internal  protrusion  whose  location  is  known.  For 
these  times,  the  nozzle  exit  stagnation  pressure  can  be  determined 
and  this,  multiplied  by  the  exposed  area  on  which  it  acted  would  give 
the  force.  The  time  at  which  the  missile  exits  the  tube,  thus  giving 
exhaust  impingement  forces  on  the  tube  face,  can  also  be  found. 

A  good  estimate  of  the  recoil-time  history  can  be  obtained  by 
plotting  all  the  separate  force-time  histories  on  graph  paper  and 
summing  to  obtain  the  total  force-time  history.  An  example  of  the 
result  of  this  process  is  shown  in  Figure  3-2  for  a  stepped  telescoping 
tube. 
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Section  IV.  LAUNCHER  SIMULATION 

4. 1  Introduction 

A  computer  simulation  program  has  been  written  which  simulates 
launcher  performance  in  the  pitch  plane.  The  program  is  capable  of 
simulating  both  smooth-bore  and  stepped  telescoping  launchers.  It 
has  provisions  for  utilizing  a  recoil  force-time  history  as  an  input 
and  for  studying  the  effects  of  thrust  misalignment  and  center  of  gravity 
offset.  Figure  4-1  is  a  representation  of  a  missile  and  launcher 
showing  the  forces  accounted  for  in  the  simulation  program  and  the 
dimensional  data  required. 

The  basic  assumption  for  the  simulation  model  is  that  the  gross 
motions  of  the  system  are  due  to  the  nonrigid  support  characteristics 
of  a  man.  The  man  has  been  replaced,  therefore,  by  a  spring  and  a 
damper.  Four  firing  positions  were  chosen  as  being  readily  simulated 
for  various  configurations.  These  positions  are  depicted  in  Figure 
4-2.  In  each  position  a  point  which  appears  to  have  the  least 
translational  motion  has  been  chosen  as  a  virtual  trunnion.  The 
spring  and  damper  are  then  converted  to  torsional  equivalents  about 
this  point. 

This  simulation  model  is  intended  only  to  provide  gross  motions. 

A  representative  example  of  outputs  from  the  program  is  shown  in  Figures 
473,  4-4,  4-5.  Results  obtained  by  the  use  of  this  program  are  predictions 
only  and  have  not  been  verified  by  test  data.  No  plans  have  been  made 
for  verification  tests  at  the  present  time. 
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:  Force  values  associated  with  blast  force  buildup. 

:  Integration  time  increment  (sec). 

:  Angular  location  of  thrust  misalignment  vector  from  zero 
reference  (counterclockwise  positive). 

:  Torque  necessary  to  induce  spin  on  the  missile. 

:  Weight  of  missile  (lbs). 

:  Weight  of  launcher  (lbs). 

:  Launcher  elevation  angle. 

:  Distance  along  X-Axis  from  missile  CG  to  front  shoe 
(Application  of  Fi)  (inches). 

:  Distance  along  X-Axis  from  missile  CG  to  rear  shoe 
(Application  of  F2)  (inches). 

:  Missile  CG  offset  along  Y-Axis  (inches). 

:  Missile  thrust  misalignment  (inches). 

:  Distance  along  Y-Axis  from  missile  centerline  to  outer 
edge  of  shoe  (inches). 

:  Distance  along  Y-Axis  from  launcher  tube  centerline  to 
trunnion  (inches). 

*  Distance  along  X-Axis  from  launcher  trunnion  to 
Application  of  Fl  (inches)* 

:  Distance  along  X-Axis  from  launcher  trunnion  to 
Application  of  F2  (inches) 

:  Distance  along  X-Axis  from  launcher  trunnion  to  launcher 
CG  (inches). 

:  Distance  along  Y-Axis  from  launcher  trunnion  to  launcher 
CG  (inches). 

:  Overall  launch  tube  length  for  smooth  bore  launcher  concept 
guidance  length  (Inside  tube  length)  for  telescoping 
launcher  concept  (inches). 

:  Outer  tube  length  for  telencoping  launcher  concept  (inches) 
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Figure  4-2.  Firing  Position  Models  for  Launch  Simulation 
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Figure  4-4.  Missile  Fitch  Rate  Versus  Time 


STEPPED  TELESCOPING  TUBE,  HOT  ROUND 
THRUST  MISALIGNMENT  INDICATED  ON  PLOT 


Figure  4-5.  Missile  Y-Velocity  Versus  Time 


3  Launch  Simulation  Equations 
Missile  Equations: 

£Me4=it.e 

I,©  -  +  Fj.E>*  +  (p3+^OC^S"kj^  -Ps(l>s-^3^ 

£F*--M,x 

M,  X  -  (Th  +Fs -P3-P4 -Vi,  Si4£A)cose  +(p,  +  rx-VJ,  cosEa)  SiO© 

M,r 

M,r  =  (fr,+ f1-uj1  cosea) cose  -  (th*fs-f,-f4-\a),  sitfEA^snj© 

Laus/cher.  Equations: 

2Mt  =  \y 

lxY=  s(^-y)  -  +F,  (i>7+x)  -Fx  (Dg--x)  +  (s+F0  (*>«.-£>*) 

-Fs(Dw-I>s)  -pct>v  +  Wi  cosCea-^)D  -Uij.siij(EA'y)Dt0 
RsLAtlVe  Motion: 

Yo-(T>w-t>jXJ“Cos  y)  i-CDa -x)  sitf  y  -CI>s-1>3Xhcoso)-I>xsin© 

r=  -y(  d  t-fcs^sioy  +  cosy-xs>cosy-(xsiNiy-©(i>y>-I>3)5M€)-ei>lCGsa 

Y  =  -  (  Du-lJs')  ¥  cos  y  -  (Dw-r>5)  y  si»oy  -J>g  tfSiN  y  +  D8$  cosy 
■*-xaf4’s»» iy  -x^cosy-zxycosy-xsiidy-CD^-D^e'coso 
-  (Dg-Da)  €►  S<£©  4- TX^NitJ©  -T>ve  obs© 
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Foe  THE  M'SStue  : 

L&t  P3=  F,  *fc 
F^.  -  Fj_*  Fc, 

TKEis/ 

I,©  =  -R£>,  *^bL+Pc-(':l+ri)Cx>s-b1)+T»(.i>+*?>3)-FS:(l)s-*3) 

FCOM  ^MALL.  ASSUMPTioN  COS©  I 

SiisJe  a  o 


Tweeepoee 

M,*-  Tw-»-Fs -UJ,*Siisl6A -F3-F4 

=  Th +Fs-V0,  S(Ki©A  -  Fc(f,+F^ 

M,Y-  -  VJt  cos  EA  +  F,  +  FX~  TH  -sin© 

I(S‘  F,C-I>,tFc(t>^-b3)]+  FXt>xtFC(l>s-l>3>ri  -tTH(tVh])3)-T^Ds-&,  ) 

L&t  e,  =  D,i-r-c(Ds-i>3) 

Bj.=  'D-X-+-  ^ 

Bj-  1>4  -«-X>3 

£4-  1>S-Dj 

T*e»a 

1,©  - -E.F.-h  +Ej7H  -  E^Fg- 
A*J|> 

F,  M  ,Y  ■♦TrtSiue  +  Vfci,  COS  6A 
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Poe  TH£  Lau^CHER  : 

-  C.'k  ■+  P|  t>-7  -r-i  be  •+  FO  X  +■  P3  4>i  +  P+4>Z 

-Py  (b«.-bS)  -Fub<.  +  ^.cos(6A-!^b9(--\Ail.Si^(£A-^b(o 

uiweee  ip  F,  is  +  >  <f>,  =  b<.-bs- 

F,  «*  -  >  4>,  =  -Cbt-Dy^ 

Ft  t*  +-  ,  4v -  bt-Dg- 
Ft  is  -  j  4 ».  r  -(Dc-bs^ 

+  F,b7  -F,.bR+tFf+Fi.')x+-FlFc4>,  +  Pi.Pc4>J. 

~{s  Clx.-t)ff)-F(.I><.+  ^cosCeA-^t^-^s.^^fcA-  *>3>,. 

I-i^-  sC'S*-^-C$-+Cb7  +  Fc4l  -+xy,  -  tb»-FC4>x-^Ft 

-Ps  Cbu-bs-VPubv-t^cosCeA-y)!^-^^  (6A-5p)b,o 

LjST  X  F**4>.  +* 

X,.-  bg  -Fc4j,-X 

i ,  -  -Fg  ( -Dc- bs)  -  FL  be  +  vV**  (sa  -  tf)  \Axsiii  (&  A  -  ^)'blo 

TheCJSFoCjE 

*  sC^-’ifJ  -c-a?  +*,f«  -  -xFt  +  V 
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FfcDK  6Q.uftTiOK)S  For.  Rei_*Tivi£  MoTic*J  ,  6eouPifJ<*  Tssms  AisJO 
L6TT»csJ4 

t-j.=  -  CE>c-r>«?)  -b®tf1'sta9  +  X^Sirttf-lXifCoS^ 

-  CBS-D3)©lc£>s©  +  'DxezsuJe 
X3  -  T>b  cos  y>  -  CDc.  -!%■  iki cos  <$ 

fc+-  Ct>s"I>3  )  Slfj© +^005© 

We  U/u/c 

Y-  *  3  &  -  X,.©  -  X  sttJ  +2rj. 

S<rjce 

F,+FLr  M,y  +TH*nJ&+  VS, cosea 

=  M ,  [ X3  -  X*e  -  X  y> +iL|  +TH  S<  n/©  -fu) , Cos  6a 

L&TT/iJ^ 

£3  =  -M ,  X.  ‘Sid  +  M  ,1^+Tn  -s inI©  +  \A cos  Gft 
S'v/e-s 

F.  +  F^r  M,X3  y  -M,X4a  +  43 


Foe  Ft  : 

i,a  - -e.f, +e2Fl  +  E3TH-e4Fs- 

E,  F,  -EiF^-  ~I,<9  +E^Th  — E+Fs1 

£*p.  +  £i*V  eiM,Xj  + 

(e,  +6,.) F,  -  £-2  M,X3  £  -  (EZMIX+  +  !,)©  +  e^+ejlH-e^Fs- 


EtM.  ..  retM.  I,  -]  6^  6j  6^ 

P| s  6  ,-t6t  ^3tf>~|_e1+Gx  ‘tel+£tJa-+6(+e1*3+‘6l+£l'FH  -  £j7glFjr 

LgtthsK 


_  _  £xM, 

£^  ~  e,+et 

i, 

e<-  =  £,+  £,. 


Ei. 

E7  -  £,+£1. 

E0=  £,+er 


£* 


P,  a  &s*3  i?  -(Esx 4+&<.)»  4  67-fc5i-e8TH-E,F£ 


^oi_ono<i  Foe  Ft  : 

Fu  s  ~e(o^3^+  (e/o^4+^t)^  +  £u  23-EftT«  +  £3F^ 

wueai*  £l0-  es-M, 

e.is  «-e7 

X5  r  £5X4+6^ 

X^  *  £(0X4 1 

^V*’’  Et2t 3  4- E3TH -G^Fg- • 

■fcg  r  £11^:3 -£sTf/ 4  £9Fy 


■iXi 
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p«  r  es  x3  ii5 

R.s -e.oX3-»  +  XLe+is 

Twee.©  acls  T*eee  casss  t©  "Be  coJsiDeaej) 

Casc  X  :  Missive  TuB6  (  tf  — 

C*%e  IT  :  Feois/T  ^Hoes  bzc>?  aRr  Cf,-o} 

C-ass  EL:  B.&A8.  ^woes  l>eop  off  C  =-F1_-o>) 

CasgT  :  &  =  © 

Iv  tf  -tfVc$  +  X,  [6sx3iy  -Xjd+iJ^QE^atf  +  Xc©  +is]+?:, 

s  r^-y)  -  /(ll.-EsX,X3-i-XfXs-6loXvX3+XtXt] 

V  -  CXs-X^y-X'inl'S+ti 
P|  =  £  s  X  j  ^  -  X  5  &  +•  2r,^ 

=  -  e,0  x3  ip  +  x*  ©  •-»■*«■, 


Cass  IE : 


F,=o 
I,  ©  =  £2_Fi  +  £r3"n4 
Fj_=  M,  X3  M, 

•*.  1,©=  -E-2M,X4&  +£tft3+e3TH 

€>  Cl,+  X4)  r  E^M,  *3^  +  6j  i-j+GjTH 

Lb  T  *9  = 

^<0  ~  M  ■•<£3 
fc7  =  £2-*3  ^E3TH 

%'S  &  -  X-IO  &  ■** 

lx  tf  -  sC*p^-y)- 

-s(^-y) -cj  +'ir(-)f1.M,X3y+Xt.M|X4d-X^3 

*-&T  X7  -  I,.  *A,.Xj  M, 

Xg  r  M,  A:t  X4 
-  ^ » “  X  t  fc  3 

•'•  X7  ip  -  sOft-ip)-  Ctf+x8e+2-t 

Souv(n><;  Poe  &  A(0D  .* 

a  =  [s(.^-ij)x,e-cxloip  +  x(eic4-x,*1]  /(k7x.9  -x8x 

-[^-y)x9-CX9a>+Xt2:7  +  ^3i(.]/0fc7>:s“^R  Jc 
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Ft=  M  X4e  +  l3 

Y  =  *3  ip  —  x4©  -x.'  suJ 

C-ASfc  IE  :  P,  - Fj_ — o 

I,®  -E3Th 
&-  I,  Th 

Iutf  +•*. 

S=ifcCs^*-^  -  c*'+*.] 

Y -X3ip-X4® -X  Sifj  +^1. 


This  launcher  simulation  program  has  been  written  for  the 
Hewlett-Packard  9100B  with  Extended  Memory  and  Card  Reader  Units. 

As  such,  it  is  composed  of  a  number  of  programs  which  perform  specific 
calculations.  Since  the  calculations  involved  were  so  lengthy  the 
basic  program  had  to  be  broken  into  two  different  parts. 

Part  I  consists  of  programs  P0,  Pi,  P2  and  ^3*  These  are  the 
input  and  initialization  parts  of  the  program.  The  input  data 
controlled  by  these  programs  is  entered  into  the  machine  via  marked 
Data  Cards.  After  data  is  input  and  initial  conditions  are  set, 

Part  II  of  the  program  is  entered.  Part  II  consists  of  programs 

p20,  P25>  p26,  p30,  p35,  p36,  p37,  p38,  p39,  p40,  p41,  p42,  p10»  P11 
and  Pi2*'  These  programs  perform  the  integrations  cf  the  equations  and 
calculate  all  output  data.  Program  P20  calculates  thrust  values; 

P25  calculates  missile  roll;  P26  calculates  motions  along  X-Axis; 

P30  calculates  friction  blast;  P35,  P36,  P37,  ?28  p39  calculate 

values  for  the  Xj  -Xio,  -  Z 5,  and  Ei  -  E31  terms;  P^q,  P41  aPd 
P42  Case  I  programs  solve  the  equations  for  the  Case  I  conditions, 

P40»  p41,  and  P42  Case  II  programs  handle  the  Case  II  equations  and 
P40>  p41,  and  p42  Case  III  programs  handle  the  Case  III  equations; 

PlO  is  the  plot  program;  Pjj  the  time  increment  program  and  Pj^  the  time 
comparison  program. 

When  the  Case  II  conditions  are  reached  the  program  will  halt 
execution  and  the  programs  for  this  condition  must  be  entered.  The 
same  procedure  occurs  for  Case  III  conditions. 

For  plots,  when  entering  Pjq  follow  this  procedure:  (1)  Enter 
PlO  in  9100  B  Memory  (2)  Press  STOP,  END  (3)  Place  switch  on 
PROGRAM  mode  (4)  Press  STEP  PROGRAM  (5)  Enter  three-digit  memory 
location  of  value  to  be  plotted  versus  time  and  (6)  Enter  Piq  in 
extended  memory.  This  must  be  performed  each  time  a  different  plot 
variable  is  desired. 
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Section  V.  CONCLUSIONS 


The  programs  presented  herein  can  be  utilized  to  obtain  an  idea 
of  the  stagnation  temperatures  and  pressures  that  may  be  imposed 
on  a  launch  tube  and  the  surrounding  area  by  a  given  rocket,  the 
internal  pressure  in  the  launcher  at  the  rc'ket  nozzle  exit,  and  the 
effects  on  missile  and  launcher  actions  of  various  dimensional  offsets 
and  misalignmerits  and  different  firing  positions.  These  programs 
are  intended  for  use  in  design  efforts  and  were  written  for  the 
Hewlett-Packard  9100B  to  enable  in-house  efforts  of  this  nature  to 
be  more  efficiently  realized. 

It  must  again  be  stated  that  the  results  of  these  programs 
have  not  as  yet  been  compared  with  actual  data  to  determine  accuracy. 
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